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ABSTRACT 

In prostate cancer, androgen receptor (AR) binding 
and androgen-responsive gene expression are 
defined by hormone-independent binding patterns of 
the pioneer factors FoxA1 and GATA2. Insufficient 
evidence of the mechanisms by which GATA2 con- 
tributes to this process precludes complete under- 
standing of a key determinant of tissue-specific AR 
activity. Our observations suggest that GATA2 facili- 
tates androgen-responsive gene expression by 
three distinct modes of action. By occupying novel 
binding sites within the AR gene locus, GATA2 posi- 
tively regulates AR expression before and after 
androgen stimulation. Additionally, GATA2 engages 
AR target gene enhancers prior to hormone stimu- 
lation, producing an active and accessible chroma- 
tin environment via recruitment of the histone 
acetyltransferase p300. Finally, GATA2 functions in 
establishing and/or sustaining basal locus looping 
by recruiting the Mediator subunit MED1 in the 
absence of androgen. These mechanisms may con- 
tribute to the generally positive role of GATA2 in 
defining AR genome-wide binding patterns that 
determine androgen-responsive gene expression 
profiles. We also find that GATA2 and FoxA1 



exhibit both independent and codependent 
co-occupancy of AR target gene enhancers. 
Identifying these determinants of AR transcriptional 
activity may provide a foundation for the develop- 
ment of future prostate cancer therapeutics that 
target pioneer factor function. 

INTRODUCTION 

Androgen signaling mediates diverse and complex func- 
tions throughout the body ranging from skeletal develop- 
ment and maintenance to spermatogenesis on ligand 
activation of the androgen receptor (AR) (1). A member 
of the nuclear hormone receptor family, AR mediates 
androgen-dependent gene expression following release 
from cytoplasmic heat-shock proteins, receptor phosphor- 
ylation and nuclear translocation whereupon AR 
homodimers bind to recognition sites within regulatory 
elements of target genes (2). This mechanism of 
hormone-dependent gene expression facihtates the fetal 
development of the prostate and regulates maintenance 
and normal function of the prostate secretory epithelium 
(1,3). In androgen-dependent prostate cancer (ADPC), 
androgen-stimulated AR function plays a vital role 
in the aberrant proliferation of epithehal cells, thus 
androgen deprivation therapy (ADT) precipitates 
marked disease regression (3,4). As is too often the case 
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however, reestablished AR activity leads to ADT resist- 
ance, marking the fatal progression to castration resistant 
prostate cancer (CRPC) (5-7). How tissue-, cell type-, and 
disease stage-specific patterns of androgen-stimulated 
gene expression are generated is an area of considerable 
research focus, and a complete picture of the activities 
determining this specificity is required for our understand- 
ing of androgen-driven prostate cancer progression 
and for the rational design of potent prostate cancer 
therapeutics. 

In nuclear receptor (NR)-regulated gene transcription, 
pioneer factor function has been studied extensively, 
providing information on the manner in which tissue- 
specific, hormone-responsive gene expression is controlled 
(8-10). Pioneer factors are characterized by their abihty to 
engage compact chromatin, initialize local chromatin 
decondensation and provide an environment amenable 
to the recruitment of transcription factors and activation 
of transcription (10). Their role in determining NR target 
gene expression was suggested on analysis of common 
DNA sequence motifs enriched within genome- wide NR 
binding sites identified using chromatin immunopre- 
cipitation combined with tiled oHgonucleotide micro- 
arrays (ChlP-on-chip) or high-throughput sequencing 
technology (ChlP-seq) (11-15). Two key pioneer factors 
emerged from complementary analyses in distinct 
hormone-related cancers, breast and prostate, driven by 
estrogen receptor (ER) and AR transcriptional activity, 
respectively. Significantly enriched within ER and AR 
binding sites were motifs recognized by the Forkhead 
box (FOX) and GATA pioneer transcription factor 
families (10,12,13,16,17). 

FOX proteins, specifically the FOXA subclass, have at- 
tracted considerable attention for their capacity to impact 
local chromatin architecture, their essential role in priming 
for temporal patterns of gene expression observed in early 
organogenesis and their influence over tissue-specific 
patterns of NR target gene expression, of particular 
interest in hormone-related cancers (9,10,18,19). In par- 
ticular, FoxAl, expressed alongside ER and AR in the 
developing and mature mammary and prostate ductal 
epitheha, respectively (20-23), has been shown to contrib- 
ute significantly to maintaining the oncogenic functions of 
these nuclear receptors in both treatment-sensitive and 
-resistant breast and prostate cancers (11,16). While the 
role of FoxAl in ER-positive breast cancers appears to be 
the positive regulation of ER-chromatin binding, resulting 
in canonical and noncanonical ligand-dependent ER 
target gene expression and breast cancer cell proHferation 
(11,24), a more complex relationship between AR and 
FoxAl exists in prostate cancer. Initial studies described 
an essential role of FoxAl in directing AR-chromatin 
binding for the activation of androgen-dependent gene ex- 
pression (16,25). However, a complementary analysis later 
revealed that FoxAl silencing resulted in the expected loss 
of many AR binding events as well the surprising gain of 
additional AR binding sites not observed in parental cells, 
suggesting both exclusionary and facihtative roles for 
FoxAl in directing AR binding events and determining 
androgen-dependent gene expression (26,27). 



While GATA family proteins have long since been 
identified as putative pioneer factors acting in concert 
with FOX proteins in establishing temporal gene expres- 
sion patterns required for fetal liver development (28,29), 
progress in understanding GATA pioneer function in NR- 
driven cancers has been largely limited to ER activity. 
GATA3 expression patterns positively correlate with 
both ER and FoxAl in the developing, mature and can- 
cerous mammary epithelium, and are in fact necessary 
for determining hormone-responsive gene expression and 
breast cancer cell proHferation (30-34). As stated, prostate 
cancer- specific AR binding sites are also significantly 
enriched in GATA motifs, and GATA2, overexpressed 
in high-risk prostate cancer, has been shown to play an 
essential role in productive AR-chromatin binding result- 
ing in androgen-responsive gene expression (12,16,35). 
Given existing evidence for the central role of GATA 
family members in defining transcription factor activity, 
the demonstrated function of GATA factors in facilitating 
hormone-dependent gene expression, and the dearth of 
information on the mechanism by which GATA proteins 
contribute functionally to NR-transcriptional regulation, 
further analyses are needed to elucidate the means by 
which GATA pioneer factors function in hormone- 
related cancers. 

In the present work, we focus specifically on the role of 
GATA2 in recruiting AR to distal enhancers elements 
of androgen-responsive AR target genes in androgen- 
responsive prostate cancer cell lines. We provide 
evidence that GATA2 acts at multiple levels to contribute 
positively to AR transcriptional activity by enhancing AR 
expression itself, facihtating AR-enhancer binding by 
establishing an accessible local chromatin environment, 
and enhancing AR target gene expression through in- 
volvement in the formation and maintenance of regulatory 
chromatin loops between AR-bound distal enhancers 
and AR target gene promoters. These mechanisms may 
account for the generally positive role of GATA2 in 
defining global AR binding that regulates AR target 
gene expression. Additionally, we reveal a complex rela- 
tionship between FoxAl and GATA2 in mediating AR 
expression and in site-specific AR recruitment facilitating 
AR target gene expression. 

MATERIALS AND METHODS 

Cell culture 

The prostate cancer cell line LNCaP was purchased from 
the American Type Culture Collection (ATCC), and 
C4-2B cells were purchased from ViroMed Laboratories. 
These cells were maintained in RPMI 1640 media 
(Invitrogen) supplemented with 10% FBS at 37°C in 5% 
C02. For individual experiments, the medium was 
replaced by phenol red-free RPMI 1640 medium contain- 
ing 5% charcoal-stripped FBS. The cells were passaged in 
our laboratory for less than 6 months after resuscitation. 

RNA interference 

Control siRNA (siControl) and siRNAs targeting AR, 
GATA2, FoxAl (ON TARGET plus^^ siRNA) were 
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purchased from Dharmacon (Dharmacon, Lafayette, 
CO). siRNA transfections were performed using 
Lipofectamine 2000 (Invitrogen, Carlsbad, CA). Twenty- 
four hours before transfection, LNCaP cells are seeded at 
a density of 5 x 10^ cells/well for six well plates or 8 x 10^ 
cells/150mm dish. The cells are transfected with 40 nM 
siRNA as described in the manufacturer's protocol and 
maintained for 72 h in hormone-free medium. Then the 
cells are harvested or treated with R1881 for an additional 
4 h before harvest. The sequences for siRNAs are Hsted in 
Supplementary Table SI. 

Western blots 

Western blot analyses were carried out as previously 
described (36). Briefly, LNCaP cells were collected and 
lysed in RIPA lysis buffer [1% NP-40, 0.1% SDS, 
50 mM Tris-HCl pH 7.4, 150mM NaCl, 0.5% sodium 
deoxycholate, 1 mM EDTA, proteinase inhibitor cocktail 
(Roch)] for 20 min on ice, and the proteins were resolved 
on 10% SDS-polyacrylamide gels before being transferred 
onto nitrocellulose membrane (Bio-Rad). The membrane 
was blocked with 5% milk powder (Bio-Rad) in 1 x TBS 
containing 0.5% Tween-20 for 1 h, washed with TBS/ 
Tween, and incubated with the specific antibodies for 
2h. Antibodies used are Hsted in Supplementary Table S2. 

Quantitative real-time RT-PCR 

Real-time reverse transcriptase PCR (qRT-PCR) was 
carried out as previously described (36). Total RNA was 
isolated from LNCaP cells with the RNeasy Mini kit 
(Qiagen, 74104). QRT-PCR was conducted on 2|ig of 
RNA by using the MultiScribe Reverse Transcriptase and 
Power SYBR Green PCR Master Mix reagents (Applied 
Biosystems), according to the manufacturer's instructions. 
Each reaction was performed in triplicate. The primers 
used are listed in Supplementary Table SI. 

Chip and ChIP sequencing 

Chip was performed as previously described (37) with a 
few modifications. Briefly, LNCaP ceUs (1 x 10 cells/1 50- 
mm dish) were plated and grown in phenol red-free RPMI 
1640 supplemented with 5% charcoal-stripped FBS for 3 
days (CeUs were treated with 1 nM R1881 for an add- 
itional 4h as required). Cells were cross-linked with 1% 
formaldehyde for 10 min at room temperature. Chromatin 
was then sonicated, diluted, and immunoprecipitated with 
specific antibodies (Supplementary Table S2) at 4°C over- 
night. Protein A-Sepharose beads were added and 
incubated for 1 h with rotation. The beads were washed 
sequentially for 10 min each in TSE I (0.1% SDS, 1% 
Triton X-100, 2mM EDTA, 20 mM Tris-HCl, pH 8.1, 
150mM NaCl), TSE II (0.1% SDS, 1% Triton X-100, 
2mM EDTA, 20 mM Tris-HCl, pH 8.1, 500 mM NaCl) 
and buffer III (0.25 M LiCl, 1% NP-40, 1% deoxycholate, 
ImM EDTA, lOmM Tris-HCl, pH 8.1) and finally two 
times with TE buffer. After being drained with a 27 G x § 
in. needle (Becton Dickinson NJ), the beads are extracted 
two times with 60 jil 1% SDS, 0.1 M NaHC03 by vortex- 
ing in a Thermomixer (Eppendorf) and pooled eluates 
were heated at 65°C for 16h to reverse the cross-linking. 



DNA fragments were purified with the QIAquick PCR 
purification kit (Qiagen 28104) and used as the template 
in quantitative PCR reactions. Power SYBR Green PCR 
Master Mix (Applied Biosystems, Foster City, CA) 
was used, and the samples were amplified with the 
StepOnePlus Real Time PCR System (Applied 
Biosystems). The primers used in real-time PCR are 
Hsted in Supplementary Table SI. ChlP-seq Hbrary prep- 
aration was performed using the Truseq ChIP sample 
preparation kit (part# 15023092) with lOng of purified 
ChlP-DNA. The libraries were amplified using 15 PCR 
cycles, and sequencing was performed on the Illumina 
HiSeq 2500 platform at the Ohio State University 
Comprehensive Cancer Center (OSUCCC) sequencing 
core (50 bp read length, single-end), with four multiplexed 
samples per lane. ChlP-seq reads were aUgned against the 
standard hgl9 build of the human reference genome with 
Bowtie vl.0.0 (38), allowing at most two mismatches. 
Reads with more than one valid alignment were 
removed, leaving only uniquely mapped reads for 
further analysis. Peaks were called using MACS vl.4.2 
(39) with P- value threshold l.Oe-8. Heatmap counts were 
generated using HOMER v4.3 (40) with bin size 50 bp. 
The raw ChlP-seq data have been submitted to the Gene 
Expression Omnibus (GEO) repository under the acces- 
sion number GSE52725. 

Formaldehyde-assisted isolation of regulatory elements 

Formaldehyde-assisted isolation of regulatory elements 
(FAIRE) was carried out as previously described (41). 
Briefly, formaldehyde was added directly to the cell 
culture medium at room temperature to a final concentra- 
tion of 1% and incubated for 10 min. IM glycine was 
added to a final concentration of 125 mM for 5 min at 
room temperature to quench the formaldehyde. Cells 
were rinsed with phosphate buffered saline and collected 
in 1.5-ml tubes. After brief centrifugation, the cell pellets 
were resuspended in 400 |il of lysis buffer (2% Triton X- 
100, 1% SDS, lOOmM NaCl, lOmM Tris-Cl at pH 8.0, 
1 mM EDTA) and incubated on ice for 5 min. Samples 
were then sonicated for 7 min (30 s on/off cycles) using a 
Bioruptor (Diagenode) at the highest intensity. The cell 
lysates were centrifuged at 13 000rpm for 10 min at 4°C 
to precipitate cellular debris, and the soluble chromatin 
was transferred to a new tube. After two rounds of 
phenol-chloroform extraction, the aqueous phase was 
combined and incubated at 65°C overnight to reverse 
cross-linking. DNA was purified with the Qiagen PCR 
purification Kit (Qiagen 28104) and eluted in 100 |il TE 
buffer. Relative enrichment in the FAIRE-treated DNA 
was calculated with DNA from untreated cells serving as 
the control. All primer sequences are listed in 
Supplementary Table SI. 

Quantitative chromosome conformation capture assay 

Quantitative chromosome conformation capture (3C) 
qPCR assays were performed as described (16) with 
some modifications. Briefly, nuclei were first cross-linked 
with formaldehyde at a final concentration of 1% for 
10 min, then digested with 400 units of Bgl II (NEB) and 
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ligated under extra-diluted conditions. After reversing the 
cross-linking, DNA was purified by phenol-chloroform 
extraction followed by ethanol precipitation. Real time 
PCR was performed using TaqMan® Universal PCR 
Master Mix (Applied Biosy stems). The probe and 
primer sequences are Hsted in Supplementary Table SI. 
The data were normalized for primer efficiency differences 
using BAC RP11-197J19 that covers the ABCC4 locus. 
GADPH loading control was used to normalize DNA 
concentration. The interaction of two Bgl II sites in the 
GADPH locus was used for comparison between different 
3C assays. 

RESULTS 

Genome-wide overlap of AR/GATA2/FoxAl binding 

To assess the global overlap of pioneer factor chromatin 
occupancy with hormone-stimulated AR-DNA binding, 
we began by performing ChlP-seq in the androgen- 
responsive prostate cancer cell Hne LNCaP. GATA2 and 
FoxAl ChlP-seq assays were performed in the absence of 
hormone, while AR ChlP-seq was performed in the 
presence of hormone to provide a picture of the extent 
to which basal pioneer factor occupancy determines 
androgen-dependent AR binding. The four lowermost 
segments of Figure lA display the density of aUgned 
ChlP-seq reads for each factor within a 10-kb window 
centered over AR binding site locations. Notably, only 
45% of all AR binding events occur in regions without 
prior pioneer factor occupancy, indicating that while add- 
itional factors Hkely facilitate AR genome-wide binding 
patterns, in this prostate cancer cell context, coordinate 
binding of GATA2, FoxAl or both is a key feature of 
AR chromatin occupancy. That combinatorial pioneer 
factor activity is a critical determinant of AR distribution 
is suggested by the observation that 13 457 sites exhibit 
AR/GATA2/FoxAl binding, while 15 342 and only 
5036 AR binding sites overlap only with FoxAl or 
GATA2, respectively (Figure lA and B). Based on these 
findings, we went on to characterize potential mechanisms 
by which these pioneer factors, in particular the lesser- 
characterized GATA2, facilitate AR binding to DNA in 
advance of androgen-responsive gene expression. 

GATA2 promotes androgen receptor expression 

A previous study reported putative GATA2 binding sites 
within the AR gene promoter in LNCaP cells, suggesting 
a potential influence of GATA2 over AR-target gene 
expression via direct regulation of AR itself (35). 
Surprisingly, our GATA2 ChlP-seq and standard ChIP 
analyses in LNCaP cells did not reveal a significant en- 
richment of GATA2 binding within this previously 
described region, constituting the promoter of the minor 
AR(A) isoform (42) (Figures 2A and B). In contrast, two 
significant GATA2 binding peaks were observed within 
regulatory elements of the full-length AR(B) isoform 
(referred to henceforth simply as AR): one 4.6 kb 
upstream of the TSS, overlapping a FoxAl binding site 
(Site 1), and the second within the AR promoter region 
(+600 bp. Site 2) (Figure 2A). While GATA2 binding to 



the promoter of the minor AR(A) isoform was previously 
shown to depend on androgen stimulation (35), ChIP 
assays for GATA2 and FoxAl in LNCaP cells revealed 
no such reliance on hormone stimulation, revealing that 
these pioneer factors engage AR regulatory elements prior 
to androgen induction of AR expression (Figure 2B). 

To evaluate the relative contribution of each pioneer 
factor to the regulation of AR expression, we individually 
reduced GATA2 and FoxAl expression by siRNA trans- 
fection and measured AR mRNA and protein levels 
before and after stimulation with the synthetic androgen 
R1881. siGATA2 transfection resulted in a significant 
reduction in AR mRNA levels prior to and following 
hormone stimulation (Figure 2C), concordant with 
the observed reduction in GATA2 recruitment to the 
AR promoter and upstream regulatory elements 
(Supplementary Figure SI A). Surprisingly, in hormone- 
depleted conditions, siFoxAl transfection resulted in 
an increase in AR mRNA (Figure 2C) despite loss of 
FoxAl binding to the upstream regulatory element 
(Supplementary Figure SIB). These results were reflected 
in the depletion of AR protein levels following siGATA2 
transfection in the presence and absence of hormone and 
in the moderate accumulation of AR protein following 
siFoxAl transfection in the absence of hormone 
(Figure 2D). Notably, while siFoxAl transfection 
resulted in an accumulation of GATA2 protein in the 
absence of hormone, this observation does not initially 
explain the concomitant increase in AR protein level, as 
no significant increase in GATA2 binding to the upstream 
AR regulatory element was observed under these condi- 
tions (Supplementary Figure SIC). 

A previous report from Cai et al. found that FoxAl- 
mediated binding of AR to a downstream repressor region 
establishes a negative feedback loop to control AR expres- 
sion (43). We performed AR and FoxAl ChIP using 
primers for the reported repressor element, and found 
that, consistent with previous results, hormone stimula- 
tion enhances occupancy of both factors to this re- 
gion, enhancing the negative feedback loop. Following 
siFoxAl transfection, basal AR occupancy of the repres- 
sor was significantly reduced in accordance with enhanced 
AR expression due to feedback loss. Hormone stimulation 
following FoxAl knockdown was sufficient to enhance 
AR repressor binding and reestabUsh appropriate 
feedback signaling (Supplementary Figure SID), explain- 
ing our observed expression patterns following pioneer 
factor knockdown. We also performed cycloheximide 
(CHX) chase assays in LNCaP cells transfected with 
Control-, GATA2-, or FoxAl -targeting siRNAs to rule 
out the impact of AR protein stabilization and/or desta- 
bilization on the observed pattern of AR accumulation, 
though only a marginal effect of pioneer factor 
silencing on AR half-life was observed (Supplementary 
Figure SIE). While additional factors must be considered 
to gain a full understanding of pioneer factor function in 
regulating AR expression, these results provide evidence 
that GATA2 is capable of binding AR regulatory elements 
independent of hormone status and is integral to both 
basal and androgen-stimulated AR expression patterns. 
In this way, GATA2 contributes positively to androgen- 
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Figure 1. Genome-wide overlap of AR, GATA2, and FoxAl . (A) Heatmap displaying the density of aligned ChlP-seq reads for AR (left), GATA2 
(center) and FoxAl (right) within a 10-kb window centered over individual and shared factor binding sites. Lettered segments indicate factor 
occupancy in four classes of AR binding sites: A— AR only, B— AR/GATA2, C— AR/FoxAl and D— AR/GATA2/FoxAl. (B) Venn diagram 
representing the proportions of shared factor binding sites relative to total factor binding sites. Lettered regions are consistent with panel (A), and 
binding site numbers are provided for each. 



responsive AR target gene expression by direct up- 
regulation of AR. 

GATA2 enhances AR target gene expression: complexity 
in combinatorial pioneer factor occupancy 

While previous studies have demonstrated the importance 
of pioneer factors in androgen-stimulated AR target gene 
expression, rigorous analyses have been limited to the role 
of FoxAl (16,26,27). As GATA motifs and indeed 
GATA2 binding events are significantly enriched along- 
side FOX motifs/FoxAl binding events overlapping AR 
occupied regions, and as GATA2 knockdown or chemical 
inhibition is understood to result in reduction of AR 
target gene expression (12,44), we chose to systematically 
evaluate the interplay between GATA2 and FoxAl in 
facilitating AR binding to androgen-responsive gene regu- 
latory elements and the impact these factors have on 
androgen-dependent gene expression. As the vast 
majority of AR/GATA2 overlapping binding sites occur 
alongside FoxAl (Figure 1), we focused our attention on 
genes apparently coregulated by these three factors. 
Referring once again to our ChlP-seq datasets, we 



identified overlapping AR, GATA2 and FoxAl binding 
events within regulatory elements of two androgen- 
responsive genes that provided informative cases of com- 
binatorial pioneer factor activity: ABCC4 and ADPGK. 

ABCC4, which encodes a member of the multidrug re- 
sistance protein (MRP) subfamily of ATP-binding cassette 
transporter proteins, has been shown to exhibit androgen- 
dependent expression (45). Importantly, three independ- 
ent studies have suggested that while ABCC4 mRNA and 
protein expression is enhanced in primary prostate cancer 
compared with nonneoplatic prostate tissue, lost expres- 
sion over time significantly correlates with metastatic 
disease progression, biochemical recurrence following 
androgen-deprivation therapy and high Gleason score 
(45-47). As such ABCC4 is a cHnically relevant 
androgen-responsive gene product with potential prognos- 
tic value in predicting the development of androgen inde- 
pendence in advanced prostate cancer. Intronic androgen- 
stimulated AR binding in the ABCC4 locus (29.9 kb 
downstream of the TSS) was found to occur alongside 
hormone-independent GATA2 and FoxAl binding sites 
(Figure 3 A and B). Basal and androgen-stimulated AR 
recruitment to this site is significantly reduced following 
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Figure 2. Novel GATA2 binding sites in the AR locus. (A) GATA2 and FoxAl binding sites in the AR locus are shown in IGV. Previously reported 
GATA2 binding site is located immediately downstream of the AR(A) TSS. (B) LNCaP cells were grown for 72 h in the absence of hormone and 
were then treated with 1 nM R1881 or ethanol for 4h before harvest. ChIP assays were performed using antibodies against GATA2 or FoxAl . 
(C) LNCaP cells were transfected with indicated siRNAs. After 72 h, cells were treated with 1 nM R1881 or ethanol for an additional 4h. Real-time 
RT-PCR was then performed using AR gene-specific primers. (D) LNCaP cells were transfected with siControl, siFoxAl or siGATA2. After 72 h, 
cells were treated with 1 nM R1881 or ethanol for an additional 4h. Western blots were performed with indicated antibodies. */'<0.05, **P<0.01 
Results of Chip and qRT-PCR assays reported as the mean of two to four rephcates, with error bars representing the standard deviation. 



siFoxAl transfection (Supplementary Figure S2A), 
and following siGATA2 transfection (Supplementary 
Figure S2B), AR recruitment is reduced nearly to the 
level observed following RNA interference of AR itself, 
demonstrating that each factor plays a central role in 
AR recruitment at this site (Figure 3C). Consonantly, 
androgen-stimulated ABCC4 mRNA expression in 
siGATA2 and siFoxAl transfected cells was reduced 
nearly to siAR levels compared with control knockdown 
(Figure 3D). To assess the hierarchical relationship 
between AR, GATA2 and FoxAl binding at this locus, 
we first characterized GATA2 and FoxAl recruitment 
to the ABCC4 regulatory element following siAR trans- 
fection. Neither FoxAl nor GATA2 enrichment was 
affected by AR knockdown in the presence or absence 
of androgen, demonstrating that both factors engage 
this site prior to androgen stimulation and independently 
of AR (Figure 3E). In addition, FoxAl recruitment was 
unaffected by siGATA2 transfection as was GATA2 re- 
cruitment following siFoxAl transfection, revealing that 
in this instance, these pioneer factors function independ- 
ently of one another in binding this regulatory element 
prior to and following hormone stimulation (Figure 3F). 

ADP-dependent glucokinase (encoded by APDGK) was 
recently reported as a novel androgen-dependent AR target 
gene with a likely role in defining cellular energy metabol- 
ism in LNCaP cells (48). While mechanistic studies of 



ADPGK have not been conducted in a model of prostate 
cancer, a general role for the enzyme in supporting glyco- 
lytic energy production under hypoxic stress conditions 
characteristic of many solid tumors has been proposed, 
and in fact, ADPGK null H460 human lung cancer cells 
exhibited decreased clonogenic survival in hypoxic condi- 
tions compared with parental cells (49). Of potential rele- 
vance to prostate cancer is the recent finding that ADPGK 
activity can induce NF-kB signaling in response to ROS 
production (50). The interplay between NF-kB and AR 
in the progression of prostate cancer has been of consider- 
able interest and one report suggests that NF-KB2/p52 up- 
regulation may be a basis for resistance to the second- 
generation antiandrogen, enzalutamide (51). Within the 
ADPGK locus, robust AR, GATA2 and FoxAl binding 
sites overlap 49.6 kb downstream of the TSS, at an 
inter genie regulatory element (Figure 4A). Here, androgen 
stimulation resulted not only in enhanced AR binding but 
also in increased FoxAl and, to a lesser extent, GATA2 
binding (Figure 4B). This is consistent with previous studies 
showing significant increases in pioneer factor binding at 
sites exhibiting adjacent AR and pioneer factor motifs 
following hormone treatment and likely reflects global 
changes in the chromatin landscape, allowing for 
enhanced factor occupancy, or stabilized pioneer factor oc- 
cupancy owing to interactions with fully assembled DNA- 
binding transcriptional complexes (12). However, this 
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Figure 3. GATA2 regulates ABCC4 expression. (A) AR, GATA2 and FoxAl binding sites are shown in IGV. (B) LNCaP cells were grown for 72 h 
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observation does raise the question of whether androgen- 
stimulated increases in FoxAl and GATA2 are a result of 
pioneer factor-like AR activity, as has been reported for the 
glucocorticoid receptor (GR), a related hormone-inducible 
transcription factor (52,53). AR knockdown revealed no 
impact on GATA2 or FoxAl occupancy in the absence 
of hormone, indicating that these factors bind independ- 
ently of AR at this location before androgen treatment. 
Yet in the presence of hormone, enhanced AR binding 
and transcription complex assembly likely stabilizes 
pioneer factor occupancy, as GATA2 and FoxAl binding 
levels following hormone treatment were diminished by 
siAR transfection (Figure 4C and D). That GATA2 and 
FoxAl act to recruit AR to this regulatory element was 
demonstrated by AR ChIP following siGATA2 or 
siFoxAl transfection, as basal and androgen-stimulated 
AR enrichment was inhibited by knockdown of these 
pioneer factors (Figure 4E and Supplementary Figures 
S3A and B). Basal and androgen-stimulated expression of 



ADPGK was also significantly inhibited by knockdown of 
AR, GATA2 or FoxAl, indicating that all three factors 
are essential for transcriptional activation at this site 
(Figure 4F). The mechanism behind this strong dependence 
was readily explained by assessing the interplay between 
GATA2 and FoxAl binding to this region. siFoxAl trans- 
fection resulted in a significant loss of GATA2 enrichment 
before and after hormone treatment, and similar results 
were obtained for FoxAl enrichment following siGATA2 
transfection (Figure 4G and H). These results indicate that 
androgen-stimulated AR binding to this regulatory element 
is facilitated following the codependent binding of GATA2 
and FoxAl to the site, and that loss of either of these 
pioneer factors results in lost binding of the complementary 
pioneer factor, binding of AR and expression of ADPGK in 
this context. 

To further support these findings, we performed a series 
of experiments in an additional androgen-responsive 
prostate cancer cell Hne C4-2B. Our results suggest that 
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elements of this mechanism of AR target gene 
coregulation by GATA2 and FoxAl are conserved 
across these two systems. A high degree of concordance 
was observed for patterns of hormone-independent 
pioneer factor occupancy of AR regulatory elements and 
regulation of AR expression (Supplementary Figure S4A 
and C). Similarly, hormone-independent pioneer factor 
occupancy was observed in advance of androgen- 
stimulated AR binding to the ABCC4 and ADPGK 
regulatory regions, and loss of either pioneer factor 
significantly abrogated subsequent AR binding 
(Supplementary Figure S4D and G). Despite the agree- 
ment between ADPGK expression changes following 
AR, GATA2 and FoxAl knockdown, siFoxAl transfec- 
tion of C4-2B cells actually induced ABCC4 expression 
while AR and GATA2 knockdown resulted in reduced 
ABCC4 expression (Supplementary Figure S4H and I), 
reveaUng cell Hne differences affecting context-specific 



pioneer factor function. This initial investigation of a 
small subset of androgen-responsive genes reveals two im- 
portant characteristics of pioneer factor function in AR 
transcriptional activation. First, the results support a 
model that GATA2 serves a generally positive role in 
mediating androgen-stimulated binding of AR to target 
gene regulatory elements subsequent to direct regulation 
of AR. Second, GATA2 and FoxAl exhibit two distinct 
models of combinatorial binding to AR target gene regu- 
latory elements: independent and codependent. 

GATA2 establishes accessible chromatin signature in 
regulatory elements 

A well-estabHshed role for pioneer factors is in facilitating 
local chromatin reorganization, allowing for improved ac- 
cessibility for additional DNA-binding transcription 
factors (10). We asked whether GATA2 might act in this 
way to enhance AR binding to distal gene regulatory 
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elements, thus allowing for androgen-dependent AR 
target gene expression. Our histone mark ChlP-seq data 
(Z.C. and Q.W., unpublished data) provided a picture of 
the chromatin environment within the vicinity of shared 
AR/GATA2/FoxAl binding sites within AR target gene 
loci of interest, illustrating that these sites are enriched 
with active enhancer-specific histone modifications 
(histone H3 lysine 4 mono- and dimethylation, and 
lysine 27 acetylation — H3K4mel, H3K4me2 and 
H3K27ac, respectively) and void of the canonically repres- 
sive histone mark, histone H3 lysine 27 trimethylation 
(H3K27me3) and the promoter-specific histone H3 lysine 
4 trimethylation (H3K4me3) (Figure 5 A and B) (54). 

We confirmed the ChlP-seq results for the ABCC4 and 
ADPGK regulatory elements, demonstrating that these 
sites exhibit typical enhancer characteristics. Relative to 
two control regions (GATA2 binding sites enriched in 
H3K27me3), androgen-independent levels of H3K4mel 
and H3K4me2 were significantly enriched, while 
H3K4me3 was only modestly enriched (Figure 5C). 
Additionally, H3K27ac, an indicator of active chromatin, 
was enriched within ABCC4 and ADPGK regulatory 
elements prior to androgen stimulation, compared to 
control loci (Figure 5D). Finally, we confirmed that 
these sites are depleted of the repressive H3K27me3, 
which showed robust enrichment within the control 
regions (Figure 5E). These results suggest that a relatively 
active chromatin environment can be found at the regula- 
tory elements of androgen-dependent genes prior to 
hormone stimulation. In addition to an active chromatin 
signature, the ABCC4 and ADPGK regulatory loci also 
exhibit enhanced accessibihty as measured by FAIRE 
(41), indicating relatively high dissociation of DNA from 
nucleosomes within these regions (Figure 5F). That 
androgen-independent pioneer factor occupancy of these 
loci correlates with this active and accessible chromatin 
signature supports the established role for these factors 
in priming AR target genes for transcriptional responsive- 
ness to hormone stimulation. 

While the mechanisms by which FoxAl influences chro- 
matin accessibility is understood to derive from its struc- 
tural similarity to the linker histone HI, allowing FoxAl 
to displace HI and disrupt heterochromatin (18,19), the 
impact of GATA2 on chromatin structure and histone 
modification status in AR signaling has been suggested 
through correlative studies, though not thoroughly 
described (44,55). Having observed an enrichment of the 
active histone mark H3K27ac, we asked what impact 
GATA2 and FoxAl have on the recruitment of histone 
acetyltransferases (HATs) to the ABCC4 and ADPGKXoci 
prior to hormone stimulation. We focused on p300, as this 
histone-modifying enzyme showed the most robust enrich- 
ment within these sites prior to androgen treatment 
(Supplementary Figure S5A and B) (56). 

p300 enrichment was significantly influenced by pion- 
eer factor status, as siGATA2 transfection resulted in 
decreased p300 recruitment to both gene regulatory 
elements in the absence of hormone, while siFoxAl trans- 
fection only significantly decreased p300 levels at the 
ADPGK locus (Figure 5G). Additionally, levels of 
H3K27ac and patterns of chromatin accessibihty at each 



gene locus mimicked p300 recruitment (Figure 5H and I). 
To our surprise knockdown of either pioneer factor 
resulted in a significant loss in H3K4mel/2 at both 
the ABCC4 and ADPGK loci (Supplementary 
Figures S5E-H). While a similar observation was made 
previously on FoxAl depletion (26,27), this is a novel 
observation for GATA2. Though the mechanism is not 
immediately apparent, pioneer factor occupancy likely 
affects recruitment of H3K4mel/2-specific histone 
methyltransferases and/or demethylases. Interestingly, 
trends in hormone-independent GATA2-mediated recruit- 
ment of p300 and estabhshment of active, accessible chro- 
matin closely resemble trends in androgen-stimulated AR 
binding (Figures 3C and 4E), suggesting a central role for 
GATA2 in providing an amenable environment for AR 
binding subsequent to hormone stimulation. That FoxAl- 
independent GATA2 occupancy of the ABCC4 locus is 
sufficient for the recruitment of p300 and the estabhsh- 
ment of active, accessible chromatin, suggests that this is 
an important mechanism by which GATA2 functions 
upstream of AR to mediate androgen-responsive gene 
expression. Together, these results demonstrate that 
GATA2, with variable input from FoxAl, plays a 
general role in enhancing androgen-stimulated AR- 
mediated expression of these target genes by engaging 
gene regulatory elements in the absence of hormone and 
recruiting histone-modifying enzymes that activate local 
chromatin. In this way, androgen-dependent gene loci 
are primed for expression following hormone stimulation. 

GATA2 facilitates regulatory chromatin loop formation 

A prominent feature of AR binding is the tendency 
toward occupancy outside of proximal promoter regions 
of its target genes. To facilitate transcription, long- 
distance interactions form between AR-bound distal 
enhancers and androgen-responsive gene promoters, 
requiring a host of intermediary protein-protein inter- 
actions involving the Mediator coregulatory complex 
(12,16,57,58). While it is understood that pioneer factor 
residence on regulatory DNA elements persists beyond the 
priming period and into periods of active transcription 
(10), their role in facihtating and maintaining activating 
chromatin loops is not well estabhshed. A previous study 
showed that FoxAl is partially responsible for estabhsh- 
ing regulatory chromatin loops in castration resistant 
prostate cancer cell lines (37), suggesting that this 
pioneer factor is able to recruit requisite loop-forming 
activities including the Mediator subunit MEDl. We 
asked whether GATA2 and FoxAl are required in regu- 
latory loop formation for androgen-dependent AR target 
gene expression by performing 3C assays (59) at the 
ABCC4 locus in LNCaP cells. 

We assessed the interaction between the downstream 
ABCC4 regulatory element, bound by AR, GATA2, and 
FoxAl, with the promoter region. The downstream site 
exhibited enrichment of MEDl prior to hormone treat- 
ment, reflected by the relatively high cross-hnking fre- 
quency observed in the absence of hormone. Following 
androgen treatment, MEDl occupancy was enhanced, 
congruent with increased cross-linking frequency 
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indicating stronger long-distance interaction between 
this site and the ABCC4 promoter (Figure 6A and 
Supplementary Figure S5A). Importantly, siGATA2 and 
siFoxAl transfection resulted in significant decreases in 
cross-Hnking frequency in the absence of hormone 
(Figure 6B), suggesting that these pioneer factors play 
important roles in facihtating basal chromatin loop for- 
mation, which may combine with basal active chromatin 
signatures and accessibihty to prime ABCC4 for 
androgen-stimulated expression. 

To assess the mechanism by which these pioneer factors 
contribute to chromatin loop formation and maintenance, 
we performed MEDl ChIP at the ABCC4 enhancer fol- 
lowing siControl, siGATA2 and siFoxAl transfection. 
Our results support previous findings that FoxAl func- 
tions in the recruitment of MEDl to distal regulatory 
elements (37), as FoxAl knockdown precipitated a signifi- 
cant reduction in MEDl occupancy. Importantly, 



siGATA2 transfection produced similar results, reveal- 
ing that GATA2 also functions in the basal recruitment 
of chromatin loop-forming activities (Figure 6C). As 
previous work has shown that phosphorylation of 
MEDl is required for locus looping and AR target gene 
expression (37), we performed p-MEDl ChIP to show 
that both FoxAl and GATA2 are central to the recruit- 
ment of this activated MEDl form (Figure 6D). It is 
worth recalHng that GATA2 and FoxAl bind independ- 
ently of one another at the downstream ABCC4 regula- 
tory element (Figure 3F); thus, decreases in MEDl, 
p-MEDl and cross-Hnking frequencies are the result of 
independent pioneer factor knockdown and reflect the in- 
dividual impact of each pioneer factor on loop formation 
and maintenance. MEDl and p-MEDl recruitment to the 
ADPGKXoQw^ following siFoxAl and siGATA2 transfec- 
tion mimicked results at the ABCC4 locus (Supplementary 
Figure S5C and D), though codependent binding of 
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Figure 6. 3C analysis of interaction between GATA2 binding site and ABCC4 promoter region. (A) LNCaP cells were grown for 72 h in the absence 
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FoxAl and GATA2 at this site prevents our evaluating 
the independent contribution of each pioneer factor in this 
instance. 

GATA2 supports genome-wide AR binding 

Having identified three prominent mechanisms by which 
GATA2 contributes to androgen-responsive AR target 
gene expression, which suggests a generally positive role 
for this pioneer factor in recruiting AR to target gene loci 
subsequent to activation of chromatin within regulatory 
DNA elements, we sought to evaluate the prevalence of 
this activity in determining AR binding patterns genome- 
wide. To address this question, we performed AR 
ChlP-seq following hormone stimulation in LNCaP 
cells transfected with Control- or GATA2-targeting 
siRNA molecules. 49 082 (79.5%) of the AR binding 
sites identified in siControl conditions were lost following 
GATA2 knockdown (Figure 7A). Among the 12 529 
(20.2%) AR binding sites that were maintained following 
siGATA2 transfection, a significant loss in ChlP-seq read 
density within these regions was observed. The loss and 



decrease of AR binding after GATA2 silencing presum- 
ably resulted from the combined loss of AR expression 
(Figure 2C and D) and GATA2-established active and 
accessible chromatin structures at these locations. 
Corroborating our single gene analyses, AR binding 
within the ABCC4 and ADPGK loci was severely 
reduced (Figure 7B). Additional, indirect effects on AR 
binding are also likely, as the nearly 80% loss in binding 
significantly outnumbers the direct overlap of the AR and 
GATA2 cistromes (Figure lA and B), though these 
effects cannot be distinguished from the upstream effect 
on AR expression. In stark contrast to the tremendous 
gain in AR binding sites previously reported following 
FoxAl knockdown in prostate cancer cell lines (26,27), 
only 131 AR binding sites were unique to siGATA2 con- 
ditions. These siGATA2-specific AR binding sites fall 
into two categories: (i) direct GATA2 exclusion of AR 
binding and (ii) trans-repression of AR binding from a 
distance. Notably only 19 sites exhibited a gain in AR 
binding at sites occupied by GATA2 in parental cells, 
indicating that >85% of gained AR binding events are 
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the result of relief of indirect or trans-repressive GATA2 
effects. Figure 7C gives examples of each of the three 
categories of AR binding sites. These results further sub- 
stantiate our claim that through its important roles in up- 
regulating AR expression and recruiting AR to activated 
chromatin, GATA2 serves a generally, if not exclusively, 
positive role in defining AR genome-wide binding 
patterns that determine androgen-responsive gene expres- 
sion profiles. 

DISCUSSION 

In the present study, we sought to close the gap in know- 
ledge between mechanisms by which prominent pioneer 
factors, GATA2 and FoxAl, contribute to AR-mediated 
gene expression in androgen-responsive prostate cancer 
cells. While FoxAl has been a major focus in previous 
studies of AR activity, limited evidence has been 
provided to illustrate a mode of action for GATA2, 
though its contribution to AR target gene expression has 



long been appreciated. Our initial genome- wide analysis 
provided compelling evidence that GATA2 plays a pivotal 
role alongside FoxAl in defining androgen-stimulated AR 
binding patterns and encouraged an in-depth investigation 
into the mechanisms used by GATA2 both in dictating 
AR binding and in modulating AR target gene expression 
(Figure 1). We showed that GATA2 operates first by 
directly regulating expression o{ AR itself via occupying 
novel upstream and promoter regulatory elements 
identified by ChlP-seq. While previous work from 
B — hm et al. reported direct regulation of AR by 
hormone-stimulated GATA2 binding at the promoter of 
the minor AR(A) isoform (35), we observed hormone-in- 
dependent occupancy of GATA2 within regulatory 
elements of the major AR(B) isoform. As GATA2 
knockdown resulted in significant loss of AR expression 
before and after androgen treatment, we suggest that 
GATA2 mediates both basal and hormone-stimulated 
expression of full-length AR from these binding sites 
(Figure 2). 
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Zaret and Carroll described a role for pioneer factors in 
"establishing competence for gene expression" by 
engaging gene regulatory elements in advance of transcrip- 
tion factors, thereby priming genes for responsiveness to 
appropriate stimuli (10). FoxAl facihtates NR-mediated 
gene expression by engaging compact chromatin, 
displacing HI to generate localized accessibihty, and also 
serves to recruit architectural components of transcription 
regulatory chromatin loops (18,37). These roles have been 
shown in various systems to contribute to NR target gene 
expression, but whether GATA2 similarly primes for 
hgand-inducible AR target gene expression has not been 
fully demonstrated. By investigating the contribution of 
GATA2 to androgen-responsive AR target gene expres- 
sion, we observed GATA2 binding in advance of hormone 
treatment and AR binding, allowing for androgen- 
stimulated AR recruitment to target gene regulatory loci 
and gene expression (Figures 3 and 4). Within our subset 
of AR target genes, we observe a positive impact of 
GATA2 on androgen-dependent gene expression, though 
it is reasonable to suspect that a repressive model of target 
gene coregulation may be adopted in a site-specific 
manner. Additionally, we have provided evidence that, 
like the documented abihty of FoxAl to preclude AR 
binding at certain genomic locations (26,27), a small 
class of GATA2 binding sites appear exclusionary of 
AR occupancy and thus also play a negative role in pre- 
venting AR transactivation in these regions (Figure 7). 
However, this subset is negligible in comparison with the 
FoxAl study, and only a minor fraction of these sites 
appear to be the result of direct exclusion of AR by 
GATA2 occupancy. Future analysis of the repressive 
mechanisms of GATA2 and/or FoxAl in determining 
AR-mediated gene expression profiles is therefore of 
great interest. 

Our analysis has provided some insight into the mech- 
anisms by which GATA2 contributes to AR target gene 
expression. Previous studies have suggested a notable cor- 
relation between GATA2 and p300 recruitment to AR 
target gene loci: a model corroborated by concurrent 
loss in binding of the pioneer factor and histone 
modifier following treatment with a natural compound 
(44). Additionally, it has been demonstrated in model 
systems that GATA2 DNA-binding and transactivation 
potential rely in part on acetylation by p300, providing 
evidence of their direct interaction (60). We have 
extended these findings to show that GATA2 has a 
central role in recruiting p300 to target gene loci, where- 
upon this HAT activates chromatin via acetylation of 
H3K27, generating local chromatin accessibihty prior to 
hormone stimulation. Within the ABCC4 gene locus, we 
can observe the individual contribution of pioneer factors 
to chromatin activation, and we see that GATA2 but not 
FoxAl knockdown results in lost p300 recruitment, di- 
minished H3K27ac levels and chromatin compaction 
relative to siControl transfection (Figure 5). Thus, we 
demonstrate that GATA2 establishes active and accessible 
AR target gene enhancers, facihtating AR binding and 
subsequent gene expression. 

We also show that GATA2, via recruitment of the 
Mediator coregulatory complex subunit, MEDl, can 



estabhsh and/or maintain basal regulatory chromatin 
loops between AR-bound distal enhancers and 
androgen-responsive gene promoters. Again, at the 
ABCC4 locus, we find that knockdown of either FoxAl 
or GATA2 independently results in significant loss of 
MEDl/p-MEDl occupancy, which destabilizes chromatin 
looping prior to androgen stimulation (Figure 6). This 
function has been previously demonstrated for FoxAl 
(37), yet the contribution of GATA2 to chromatin loop 
formation has been more elusive. Occupancy of the 
GATA family member GATAl along with its coactivator 
FOGl at the P-globin locus has been shown to correlate 
with loop formation (61), providing the most definitive 
evidence to suggest a direct involvement of GATA 
factors in estabhshing and/or maintaining chromatin 
loops. While additional factors may be identified that 
can stabilize MEDl -enhancer occupancy and basal chro- 
matin loop formation, we suggest that androgen-inde- 
pendent GATA2 occupancy of AR target gene 
enhancers is a critical determinant of MEDl -mediated 
locus looping in preparation for hormone-stimulated 
gene expression. Previous studies have implied significant 
roles for various transcription factors and coactivators in 
loop formation (58,61-64), and our current findings 
suggest that GATA2 plays a central role in these processes 
via recruitment of p300, allowing for basal chromatin re- 
modeling, and MEDl, facilitating basal locus looping. 

Our genome-wide analysis of the impact of GATA2 
silencing on AR distribution reflects two key functions 
of GATA2 in these prostate cancer cells. First, GATA2 
knockdown significantly reduces AR protein levels and, 
subsequently, AR ChlP-seq signal intensity across the 
genome. Second, GATA2 silencing results in the loss of 
active chromatin signatures in regions primed for 
androgen-stimulated AR binding in parental cells, thus 
global ChlP-seq enrichment is diminished. While the com- 
prehensive impact of GATA2 loss on transcription- 
regulatory chromatin looping remains in question, we 
have provided evidence that the loss of this critical 
GATA2 function, in combination with the dramatic 
impact on AR binding patterns, will broadly influence 
AR target gene expression. 

Our findings support and further contribute to the 
growing body of knowledge of GATA family pioneer 
transcription factors that facilitate nuclear receptor- 
mediated, hormone-responsive gene transcription. 
Demonstrating that GATA2 is capable of facilitating 
open and active chromatin prior to hormone stimulation 
is concordant with previous observations that nucleo- 
some-depleted regions (NDRs) in AR target gene enhan- 
cers overlap GATA2 binding sites that appear to 
be necessary for NDR maintenance in the absence of 
androgen (65). Additionally, the ER-collaborating 
pioneer factor, GATA3, important for breast cancer 
gene expression profiles, exhibits binding patterns that 
correlate with genome-wide ER-anchored chromatin 
loops associated with estrogen-responsive gene transcrip- 
tion (62,66). Further analysis is required to reveal whether 
or not GATA2 shares common functionality with 
GATA3, which appears capable of activating and/or 
silencing local chromatin in a site-specific manner 
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leading to enhanced and/or inhibited chromatin loop for- 
mation (66). 

In summary, we have provided evidence that GATA2 
contributes positively to androgen-dependent AR target 
gene expression by three distinct mechanisms. First, 
through binding regulatory elements of the AR locus, 
GATA2 directly upregulates basal and hormone 
stimulated AR expression. Secondly, by recruiting the 
histone acetyltransferase p300, GATA2 facilitates basal 
AR target gene enhancer activity and promotes chromatin 
accessibility, allowing for subsequent androgen-stimulated 
AR binding and target gene expression. Finally, GATA2 
recruits MEDl/p-MEDl to AR target gene enhancers 
prior to hormone stimulation, leading to basal chromatin 
loop formation/maintenance at androgen-responsive gene 
loci. These findings encourage future focus on elucidating 
additional means by which GATA2 contributes to AR- 
mediated gene expression to further our understanding 
of how this pioneer transcription factor influences AR 
binding patterns and downstream target gene expression 
in prostate cancer. Revealing critical determinants of 
AR activity offers an important opportunity to direct 
new therapeutic strategies at upstream components of 
the androgen signaling axis. As current prostate cancer 
therapeutic modalities too often fail in preventing bio- 
chemical recurrence in advanced cases, focusing future 
treatments towards activities of factors like GATA2, 
correlated with advanced prostate cancer (35), may 
prove a potent strategy. 
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